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Abstract 

The long standing contrast between Boltzmann’s and Gibbs’ approach 
to statistical thermodynamics has been recently rekindled by Dunkel and 
Hilbert jT|, who criticize the notion of negative absolute temperature (NAT), 
as a misleading consequence of Boltzmann’s definition of entropy. A different 
definition, due to Gibbs, has been proposed, which forbids NAT and makes the 
energy equipartition rigorous in arbitrary sized systems. The two approaches, 
however, are shown to converge to the same results in the thermodynamical 
limit. A vigorous debate followed ref. [I], with arguments against miii and in 
favor mmm of Gibbs’ entropy. In an attempt to leave the speculative level 
and give the discussion some deal of concreteness, we analyze the practical 
consequences of Gibbs’ definition in two finite-size systems: a non interacting 
gas of N atoms with two-level internal spectrum, and an Ising model of N 
interacting spins. It is shown that for certain measurable quantities, the 
difference resulting from Boltzmann’s and Gibbs’ approach vanishes as TV -1 ' 2 , 
much less rapidly than the 1/N slope expected. As shown by numerical 
estimates, this makes the experimental solution of the controversy a feasible 
task. 
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Preamble 

A first version of the present work was delivered to ArXiv months ago, and, 
jointly, to The Physical Review E. After a long process, the manuscript was finally 
rejected, due to the negative opinion of one referee (the other was, istead, very 
favorable), concerning the application of Gibb’s entropy to Weiss ferromagnetism. 

ArXiv is an open space, where people can freely show their results and opinions, 
leaving any judgement to a free audience of readers. Hence, it seems to the author 
that the above piece of information is important: the reader is now acquainted that 
there was a contrast, on the content of the work, and, after a reading, she/he can 
judge about the terms of the contrast, and decide which opinion is right (if any: 
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comments are welcome). Actually, sometimes the debate raised by a work is more 
fruitful than the work itself. 

At the end of the work, a post-amble is added, in which the reasons of the 
referee’s negative judgement are reported (as far as the author did understand 
them). 


1 Introduction 


The notion of negative absolute temperature (NAT) [8, (9] has been recently re¬ 
freshed by Braun et al uni . The main factor of novelty was the claim that the 
NAT regime was attained in a system with a continuous upperly bounded spec¬ 
trum, in contrast to the preceding experiments based on two-level systems 019]. 
Soon after, however, the notion of NAT was criticized by Dunkel and Hilbert [1], 
as a misleading consequence of Boltzmann’s definition of entropy 


S B (E) = In 


Sh M, e 
. 'n 


(i) 


in an isolated (microcanonical) system of Hamiltonian H(rf), r) being any state 
variable and 6 a generalized Kronecher symbol (Boltzmann constant ks = 1). 
Actually, from the definition of absolute temperature 


T = 



( 2 ) 


it is seen that NAT occurs whenever S(E) is a decreasing function of the energy, 
which is possible for the entropy 0, if the energy is upperly bounded (i.e. H(rf) < 
Em ax) and the higher energy levels can be overpopulated by suitable external 
processes. As stressed in Ref. 1], instead, Gibbs’ entropy 


S G (E) = In 


J2®(E-H(ri)) 

- V 


(3) 


(©(•) being the Heaviside function) is an increasing function of E and thereby ex¬ 
cludes any NAT regime by definition. In Ref. (T] the authors stress that Gibbs’ 
definition of entropy refers to a microcanonical system (a claim that will be recon¬ 
sidered in what follows) and satisfies the equipartition theorem in any case, while 
Boltzmann’s entropy fails in very small systems, with a number of particles of order 
unity. 

What we playfully call the ” Boltzmann vs Gibbs match’j^] is right the debate 
about Boltzmann’s (eq.n 0 ) and Gibbs’ (eqn. 0 ) definition of entropy mmm 

1 This title should not be confused with Jaynes’ one EH: Gibbs vs Boltzmann entropies, that 
deals with the (supposed) inadeguacy of Boltzmann combinatorial method in describing interact¬ 
ing systems. 
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G3IB10- As it is a rule for most polemics involving the entropy, the discussion 
may look a little bit academic, especially in view of equation (14) of ref. pQ, which 
shows that the two entropies ([!]) and ([3]) give the same temperature expression, 
in the thermodynamic limit, with a difference that vanishes as the inverse heat 
capacity. Given that Boltzmann and Gibbs’s picture lead to the same results in 
large systems, the alternative looks as follows: either there exists a measurable 
size effect displaying some difference between eqn.s ([l]) and ([3]), or any discussion 
remains confined to a merely speculative level. The second possibility looks, at 
a first sight, much more likely, once one efforts the question how small the size 
should be, for this effect to be detectable. Actually, one might be tempted to 
assume that the smallness criterion is determined by 1 /N (N being the number of 
particles), since from eq.n (14) of ref. [1J the temperatures resulting from eq.n ([3]) 
and ([I]) differ by terms proportional to the inverse of an extensive quantity (the heat 
capacity). If this were always the case, any experimental test, like one based on 
the ’’minimal quantum thermometer” suggested in ref. |Tj, would become extremely 
difficult, if not impossible. Fairly surprisingly, we shall show that, in certain two- 
level systems, the difference between certain measurable quantities, derived from 
eq.n (J3| and ([!]), vanishes as 1/y/N. This result makes the experimental comparison 
between the two approaches much more accessible, since detecting effects small to 
order 1/y/N is certainly easier than to order 1/N, especially in view of the recent 
technical progresses of small particle physics. This is the ’’finite-size match” that 
we are going to outline in the next sections. 

Though we shall refer to two-level systems, that represent the preferred frame¬ 
work for NAT experiments, we are not specifically interested in the NAT problem. 
We address to a more fundamental point, that intrigues people since the founda¬ 
tions of statistical thermodynamics, i.e. the definition of entropy in terms of the 
microscopic dynamical states of the system. Very recently, Hilbert, Hanggi and 
Dunkel performed an exhaustive survey of several possible definitions of entropy in 
continuous spectrum systems [12] , with special attention to their consistency with 
the three principles of thermodynamics. It turns out that the only definition of en¬ 
tropy that satisfies the three principles in arbitrary small systems, is eq.n ff The 
flaws resulting from other definitions are, in turn, finite-size effects whose measured 
presence would support a widespread opinion (contrasted by Hilbert, Hanggi and 
Dunkel) that thermodynamics applies only to large systems. 

Another point that we address in the present work is the role of the con¬ 
straints in the construction of the thermodynamical functions. Actually, differ¬ 
ent constraints (microcanonical, canonical, grandcanonical) yield different fluctu¬ 
ations about the equilibrium values of the thermodynamical quantities. Another 
widespread opinion is that the equilibrium values are, themselves, independent from 
the constraints. This justifies the experimental evidence that the constraints are 
irrelevant in the thermodynamical limit (and far from the phase transitions), since 
the relative weights of the fluctuations vanish as 1/y/N. However, the difference 
between eq.n ([3| and ([l]) is not due to the fluctuations, but to the genuine notion of 
thermal equilibrium. Indeed, Gibbs’s entropy ([3]) postulates what, in Boltzmann’s 
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languag^J sounds like a sort of new ergodic hypothesis, that (in contrast to what 
claimed in ref. [T]) looks hardly appliable to micro-canonical systems. Why should 
states with energy less than E be involved, if E is strictly conserved? In this con¬ 
cern, it is only Boltzmann’s entropy 0 that is consistent with the dynamics of an 
isolated system. A reasonable possibility is that Gibbs’ and Boltzmann’s entropies 
refer to canonical and micro-canonical constraints, respectively, which is right the 
opposite of what the supporters of Gibbs’ entropy do claim The difference 

between canonical and microcanonical systems will come into play in the study of 
Weiss ferromagnets (Section [3]) . All the way, it leaves a relevant question pending, 
to which finite-size experiments could give an answer. 

What precedes motivates our program of giving the size effects of interest an 
analytical form, then studying the feasibility of ad hoc experiments revealing which 
of the two entropy definitions does fit better with physical reality, both in canonical 
and in micro-canonical systems. In order to avoid any possible factor of confusion, 
it is worthwhile stressing the strategy used in the next sections. We calculate the 
quantities of interest as functions of the number N of particles, according to Gibbs’ 
definition eq.n ([3]). Then we compare such quantities with those obtained according 
to Boltzmann’s definition ([I]), under the same condition of large but finite N. In 
Section [3] (the Ising ferromagnet) this procedure is manifest (see, for example, Fig.s 
2, 3 and 4). In Section [2j the Boltzmann-formulated quantity (the energy of the 
two-level gas) is reported from other papers (in particular, Ref. [Ij). 

Another caveat is in order: the current way to represent the differences between 
Boltzmann’s and Gibbs’ entropies is writing the Boltzmann/ Gibbs ’’temperatures” 
Tb,g(v)^ as functions of other state parameter(s) r). However, this does not mean 
introducing different Boltzmann/Gibbs ’’thermometers”, but different relationships 
between the measured temperature T and 77 . Right to point out that T is the same 
measured temperature , both for Boltzmann and for Gibbs, we will not speak about 
Boltzmann’s or Gibbs’ ’’temperature” (apart from rare cases) and will not append 
any subscript G or B to T(rj). 


2 A gas of non interacting two-level particles 

The simplest system for testing the differences between the two entropies Q and Q 
is a non interacting gas of N two-level particles, with populations n± in each level 
±e. While Boltzmann predicts n+ = n_ = N/ 2, for T —» Gibbs predicts n + ( T) 
increasing continuously up to an inverted population regime n + —> N, n_ —> 0 at 
T = 00 . At this stage one might wonder why laser devices need a resonant external 
field, in order to get an overpopulated level, if this could be obtained simply by 
thermal activation. The point is that in the Gibbs framework, the size of the 
gas (represented by the particle number N) plays a crucial role in determining how 

2 Which is hated and rejected by many Gibbs supporters E31- 

3 We implicitly assume that the two levels are non degenerate or have the same degeneration. 
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high the temperature should be, to achieve a significant overpopulation of the upper 
level. This is what we are going to show in the present section. 

On setting n± = (TV ± m)/2, the energy of the gas (as far as the two-level part 
is concerned) reads: 

£(m) — e ui. (4) 

Since the condition £(m) < E implies m < E/e, Gibbs’ entropy reads: 


S G (E) = hr 


E 

m<E/e ^ 2 


N\ 


N-\-m ^ | ^ N—m ^ | 


If TV 1, the entropy can be expressed in an integral form, making use of Stirling 
formula for the factorials: 


S G (E) = In 


~>N 


with 

and 


J -1 

E 

N~e ^ 


= <t>(z) + In ( x ~ z ) + 1 g2 ( x - z f + ■ ■ ■ 


dxe-Z+M , 

(5a) 

[-1.1] 

(5b) 

x)ln(l — a:)] = 




(5c) 


The series expansion in eq.n (5c I will be used in Appendix A. The temperature 
follows from eq.n (|2]) , on account of eq.ns (J5|: 

r z(E) 

T{E) = Ne e^ (i(B)) J dx e~^ x) . (6) 

For TV Tg> 1, the following expressions are derived in Appendix A: 


T(E) 

2e 


In 


1 - z(E) 
1 + z{E) 


i -i 


for E < 0, \E\ » e\j — (7a) 


1 irN 

2 V 2 


for E = 0 


(7b) 


1 TV N 

r tv ,, i 

v^” exp 

-2<f>(z(E)) 


1 + erf z{E) 


for 0 < E <C e TV . 


(7c) 
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The condition E < 0 corresponds to the thermal regime in which the Boltzmann 
temperature Tg is positive and the lower level is always more populated then the 


higher one. Actually, equation (7a) recovers exactly the standard result: 


z(T) = 


e -t/T _ e e/T 
e e/T + e -e/T ’ 


( 8 ) 


as can be easily shown by solving (7a) with respect to z(E), then using eq.ns (5b) 
and Q. The differences between Boltzmann and Gibbs at finite N becomes di- 
vergingly large in the range of critical energies: 


- ey/W/2 < E < 0 . 


(9) 


A E 


Indeed, it is well known that \im E ^ 0 -T(E) = oo, according to Boltzmann, for 
any value of N (see, for instance, Fig. 1 in ref. [I]). In contrast, equation (7b) 
shows that lung^o- T{E ) = Tq oc y /TV is finite , but diverging with \/N, according 
to Gibbs. The l/-\/iV-discrepancy between Boltzmann and Gibbs is expressed by 
the ratio AE/(eN), between the range of critical energies (eq.n ([9])) and the total 
energy range. Note that the condition E > 0 corresponds to what would be the 
NAT regime, that should be attained by ’’super heating” the system above E = 0, 
according to Boltzmann’s picture^ 

From a formal viewpoint, the results eq.ns 0 , obtained in the Gibbs framework, 
cure the singularity of the Boltzmann temperature at E = 0 and bring the gas 
continuously to any positive value of the energy, i.e., to any thermally activated 
overpopulation of the upper level. However, it is easily shown from eq.n ( |7c[ ) that 
the temperature diverges exponentially with N, for E > 0 (Fig. 1), while the 
crossing temperature 


To = ^ (10) 

at E = 0 (eq.n ([7b|), diverges with the square root of the system size. Hence, the 
regime E > 0 would be unaccessible to any thermal process, for N ^ oo, even in 
the Gibbs framework. In this (non trivial) sense Gibbs and Boltzmann do converge 
to the same results in the thermodynamic limit. However, in a finite-size system, 
the regime E > 0 (or, equivalently, z > 0) is not forbidden in principle, according 
to Gibbs, and could be accessed by an appropriate experimental set up. This is 
what we are going to exploit in what follows. 

Let us derive an expression for z as a function of T and N, recalling that z 
(eq.n ( |5b[ )) is the fractional difference between the higher and lower level popula¬ 
tions. Assuming 0 < z € 1, which means <f>(z) = z 2 + - ■ ■, and defining 9 = z^/N/2, 

4 This is what Sokolov criticizes as ’’hotter than hot” Yj. 


6 











Figure 1 

ln(T/e) 



Figure 1: Gibbs temperature of a two-level gas in the Boltzmann NAT 
regime E > 0. 

The logarithmic plot shows the exponential increase of T in the particle number N, 
for different positive values of z, corresponding to a thermally activated overpopu¬ 
lation of the upper level. 
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equation (7c) yields: 


9 = 


In 


T 


T G (1 + erf(0)) 


1/2 



(T > 2 Tq) , 


( 11 ) 


since erf(0) ranges between 0 and 1 and does not affect the result significantlj0 
At fixed temperature T (> 2 Tq), the overpopulation factor z decreases with the 
square root of the system size. Due to the logarithmic dependence, the value of z is 
almost insensitive to T. Hence, the feasibility of an experiment possibly measuring 
0 > 0 at a positive temperature rests on two main difficulties: first, a crossing tem¬ 
perature Tq oc y/N not so large as to prevent high precision measurements; second, 
an overpopulation factor z oc y/l/N not so small as to escape the instrumental 
sensitivity. 

Consider a gas of N atoms of 3 He, in a volume V, each carrying a spin magnetic 
moment psohr ~ 10 -2 °emu. Let B be an external uniform magnetic field. From 
eq.n ©> the total induced moment ptot = ±PBohr {n+ — n_) reads, according to 
Gibbs: 


t^tot 


—sign (B) pso hr 



(12a) 


while Boltzmann (eq.n ([8]) with e<T) yields: 


Mtot ~ sign (B) p,BohrN^ , 


(12b) 


with e = \B\ /i Bohr• Note that p,f ot is antiparallel and pf ot parallel to B, since the 
former comes from the overpopulation of the upper level, while the latter comes 
from the overpopulation of the lower one. SQUID magnetometers, that are the 
highest sensitivity instruments currently available, can measure magnetic moment 
intensities down to fj, m « 10 _8 emu. For T larger than, but comparable to Tq, the 
condition \pf ot \ ~ (eq.n (12aI) yields the order of magnitude N « 10 24 for the 
minimum number of atoms required to produce a detectable value of p,f ot . Efficient 
SQUID devices usually operate at temperatures of 10 K at most, which determines 


an upper limit for Tq. From eq.n (10), however, this turns into an upper limit for 


\B\ = 


2 Tq 
N 7 T PBohr 


10~ 7 G. 


(13) 


With those values of N, B and T, one sees that \pf ot \ ~ 10~ 8 emu (eq.n (12bI) 
is comparable in magnitude to so that the measurable difference between 

Boltzmann’s and Gibbs’ predictions results mainly in the opposite orientation of a 
small induced magnetic moment. This effect is l/-\/]V-small because the antiparallel 


5 The condition T > 2 Tq ensures that the quantity in square root is positive, for any value of 
9, but the solution of the equation exists for any T > Tq. Note that 2<1 does not imply 6 <C 1. 



















induced magnetic moment p^ ot increases as \/N (eq.n (12a)), while the parallel 
induced moment pf ot is proportional to N (eq.n (12b])). 

As for the volume V of the gas, one should recall that the size of the sample 
does also play a role, in the feasibility of the experiment. It can be easily seen that 
V ~ 10 cm 3 , corresponding to a density of 10 23 cm -3 of 3 He atoms, at T ss 10 K, 
yields degeneration effects less than 10%. This looks a still tolerable error factor for 
the implicit approximation assumed so far, that the two-level energy is statistically 
independent from the translational energy (which is false, for degenerate gases). 

In conclusion, some cube centimeters of 3 He, with density about 10 23 cm -3 , at 
a temperature about 10 K, under the action of a magnetic field of 10~ 7 G, could 
be a good candidate as a referee of the Boltzmann vs Gibbs match. If the mea¬ 
sured magnetic moment (spin polarization) of the gas, induced by the field, were 
anti-parallel to the field itself, the resulting thermally activated overpopulation of 
the upper level would provide a strong experimental support to Gibbs’ definition 
of entropy. On suitably changing the experimental set up, the hypothesis that 
Gibbs’ and Boltzmann’s entropies do refer to different constraints (microcanonic or 
canonic) could be exploited in turn. 


3 Weiss ferromagnetism revisited 

Another system for testing the consequences of Gibb’s definition of entropy is a 
ferromagnetic material, modeled by an Ising lattice of N = n+ + n_ interacting 
magnetic moments with only ”up” (+) and ’’down” (—) orientations. In this case 
one has to face the problem of applying Gibbs’ picture to a phase transition (para¬ 
magnetic O ferromagnetic), which turns out to be a non trivial issue. 

In the mean field approximation, the energy of the Ising model reads: 

£(z, P) = ~JN (z 2 + 2pz) , z = H+ N n ~ , (14) 

where J > 0 is half the coupling constant of the moment-moment interaction, 2 
is the magnetization and p = poB/(2J) is the energy contribute from an external 
uniform magnetic field, in units of 2 J. In a first-order approximation in p (|p| <C 1), 
the condition £(x, p) < E implies: 

x > |z(.E)| — p or x < — \z(E)\ — p , 


with 


\z(E)\ = 


E 


JN 


1/2 


and - (1 + p) < ~ < 0 . 


In the same limit N ^$> 1 and with the same method used for eq.n ([3], the preceding 
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formulas yield, for Gibbs’ entropy: 



Figure 2: Gibbs’ and Boltzmann’s temperatures of a Weiss ferromagnet 
(Ising model) in zero external field (N = 200). 

Dashed and full lines refer to Boltzmann’s and Gibbs’ temperature, respectively, as 
functions of the relative magnetization intensity \z\. For the meaning of z we i SS , z* 
and zq , see the text. 


Let the external magnetic field vanish (p = 0). On applying eq.n ([2]) to eq.ns (14) 
and (15), the relationship between Gibbs’ temperature and the relative magnetiza¬ 
tion intensity \z\ is shown by the full line in Fig. 2. The dashed line, instead, refers 
to Boltzmann’s temperature. The right side of Fig. 2 shows that both T-curves 
initially increase with decreasing \z\, which is what one expects for a ferromagnet 
whose spontaneous magnetic order tends to be destroyed by an increasing thermal 
disorder. Boltzmann’s temperature increases on and tends to the Curie tempera¬ 
ture T c = 2 J for \z\ 0. Gibbs’ temperature, instead, attains a maximum Tm, 

lower than T c , then decreases with decreasing \z\, down to zero, at \z\ = 0 (left side 
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of Fig. 2). A common drawback of both curves is that they are unable to account 
for the paramagnetic phase at T > T c . In addition, Gibbs’ temperature exhibits 
further intriguing aspects: for each T < Tm, there are two possible values z we i SS 
and zq of the relative magnetization intensity. Since 8T / dz is positive in zq (and 
negative in —zg), according to equation ( fl4| ), it is easy to show that the specific 
heat at constant volume is negative in ±z<Tj 

<0. (16) 

±ZG 

Furthermore, d£/d\z\ is finite at |z| = z*, in which the derivative of T(\z\) 
vanishes. Hence: 



lim| 


|-h 


88 

W\ 



(17) 


Equation shows that the specific heat Cq diverges positively or negatively 
for T —► Tm , depending on which state, z we i SS or zq, is adopted. Though a phase 
transition may well produce a multiplicity of equilibrium states at the same T, V, ■ ■ ■, 
in this case it is clear that the ’’states” ±zq are spurious and that the range of 
values 0 < \z\ < z* in Fig. 2 is unphysical in some sense. Here we suggest an 
argument that provides a plane interpretation of Fig. 2, but challenges the position 
that Gibbs’ entropy and temperature refers to microcanonical systems only, while 
Boltzmann’s picture refers to canonical systems Q On assuming that eq.n ([2| does 
apply, at least locally, one has: 


(88 N 


'{88/dri)y 

\8S, 

K 

.( dS/8r]) v _ 


for any state parameters rj that characterizes the system, and attain the equi¬ 
librium values r] e q. The next step follows from rephrasing the preceding formula 
as: 


d(£- TS ) 


dg 


= 0, 


e,/ y T 


(18) 


which yields the extremants of the Helmholtz free energy ^ = £ — TS. In 
Appendix 2, it is shown that the spurious states zb zq correspond to maxima of 


6 In the next formulas we use the following convention: \dA/drj\ means: derivative of quantity 
A(r]), calculated in rj eq \ (dB/dX ) v ... means: derivative of quantity B with respect to X, keeping 
V, • • • constant. 

7 In ref. m, for example, it is shown that Gibbs’ entropy is the unique expression that fits 
rigorously with Clausius definition of temperature in a microcanonical system. When dealing 
with a canonical system, instead, it is Shannon’s entropy (Shannon’s and Boltzmann’s entropy 
coincide, in a canonical system) that shares the same property. In ref. [4] it is stressed that 
the difference between Sb and Sg could result right from the difference between canonical and 
microcanonical constraints (See footnote [6] in ref. El) 
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the Helmholtz free energy T q = £ — TSq, while the state z = 0 is always a 
cuspid minimum, i.e. \I>g(z) — ’I'g(O) = T\z\\J 2 / {N tt) for small \z\ (eq.n (22)). 
However, this makes sense only if the temperature T, defined through eq.n ([2j), 
refers to a canonical system, for which the minimization of the Helmholtz free 
energy is the correct approach to thermal equilibrium. Hence, the interpretation 
just outlined stems from the assumption that Gibbs’ entropy does apply to canonical 
systems too (or only), which is a pending question. Leaving the alternative between 
microcanonical or canonical use of Sg to future discussions and investigations, we 
take the states z we i SS (with positive specific heat) as the ones corresponding to the 
true thermal equilibrium, and stress the most relevant result following from Fig. 2: 


Figure 3 

Cy/N 



Figure 3: Gibbs and Boltzmann heat capacities as a function of T in a 
Weiss ferromagnet (N = 200). 

A finite discontinuity at T c = 2J (dashed line) vs a divergent discontinuity at 
Tm < T c (full line) characterize Boltzmann’s and Gibbs’ heat capacity, respectively. 


Figure 3 shows the behavior of Boltzmann’s (dashed line) and Gibbs’ (full line) 
heat capacities: while the former displays a finite discontinuity at the Curie tem¬ 
perature T c = 2 J (which is reported even in the pedagogical literature |14jl. the 
latter exhibits a divergent discontinuity at the critical temperature Tm, lower than 
T c . 

As far as our aims are concerned, the question turns back to the possible exper¬ 
imental evidence of the diverging heat capacity mentioned above. At a first sight, 
this effect looks huge and very easy to observe. Indeed, one might wonder why it 
has not been reported before, if Gibbs is right. The problem is much less trivial, 
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actually. First, it can be seen that \T C — Tm\ oc 1 /y/N, i.e. the difference between 
the standard Curie temperature and Gibbs-Curie temperature Tm is a size effect, 
vanishing as the inverse square root of the system size (see Fig. 4). This means that 
the temperature range in which the diverging behavior could be observed vanishes 
in turn, with at least the same l/v^-slope. In addition, under the hypothesis that 
Gibbs’ picture is only microcanonical, the Cb divergence could be observed only if 
the phase transition is approached in isolation. In ref. [6], the practical difficulty 
of implementing such a process is stressed with some interesting details. So, the 
experimental check of Fig. 3 turns out to be far from easy. If, instead, Gibbs’ pic¬ 
ture were also canonical, Section [4] outlines a possible experimental route, based on 
magnetic measurements, that turns out to be feasible, with appropriate high-level 
instruments. 


Figure 4 

ln(l-7Wr c ) 



ln(N) 


Figure 4: Log-Log plot of 1 — Tm/T c as a function of the system size N. 

The numerical data (full line) are compared to the predicted 1 /y/N slope (dot- 
dashed line). 


4 Weiss Ferromagnetic transition as a canonical 
process 

As anticipated in Section [3j here we assume that both Boltzmann’s and Gibbs’ 
entropies apply to canonical systems, and thereby proceed to minimize the free 
energies 


* a {z)=£(z)-TS a {z) (a = B,G ) (19) 
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with respect to the magnetization z, according to eq.n (181. For the sake 
of brevity, we do not report the explicit formulas for which is a standard 

pedagogical issue. As for \I>g( 2 ), according to eq.ns (14) and (15), one gets, to first 
order in p : 


p ) = — JN (z 2 + 2 p z) — T In 


2 N+1 / dxe“^ x) 
J\z\ 


( 20 ) 


First, let us study the perfectly symmetric system at p = 0. The extremants of 
^g{z, 0) follow from the equation 


d*g{z, 0) 
dz 


= -2 JN z + sign ( 2 ) T 


which yields (see Appendix A): 


M T , ( 1 + w 
w = 37 ln TVu 


dx 


'Id 


for I z I > 


= 0 . 


2N ’ 


(21a) 


T 


JV2N~i 


= z G 


for \z I < 


(21b) 


In agreement with the rule that Gibbs and Boltzmann converge to the same results 
in large systems, Boltzmann’s picture of Weiss ferromagnetism is recovered from 
eq.ns (21), in the limit N —> oo. Indeed, one easily sees that in this limit 2 = 


0 is always a solution, while other two solutions ±z we i SS exist, from eq.n (21a), 
provided T < T c = 2 J. The solutions ±z we i SS correspond to two minima of the 
Helmholtz free energy and yield the spontaneous magnetization, below the Curie 
temperature T c , where, according to Boltzmann, z = 0 is a maximum (i.e., an 
unstable equilibrium state). Above T c , the solution 2 = 0 is the only minimum of 
Boltzmann’s free energy, which marks the transition from the ferromagnetic to the 
paramagnetic phase. All this is displayed by the dashed line plots in Fig. 5. 

In a finite system, Gibbs and Boltzmann pictures of Weiss ferromagnetism differ 
for small values of | 2 |. Indeed, for T < T M (recall Fig. 2), the solutions ±zq = 

±T/ (^J\Z2Ntt^ of eq.n (21b) correspond to two maxima of 4 >g(z, 0), while 2 = 0 
is always a minimum, as can be seen by the full line plots in Fig. 5 and by noticing 


that (eq.n (21b)): 


lim z _> 0 ± 


d*g(z, 0) 
dz 


= ±T 


2 

Nn 


( 22 ) 


which shows that the minimum at z = 0 corresponds to the terminal point of a 
downward cuspid, as sketched in Fig. 5. 
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Figure 5 a 



Figure 5b 



Figure 5: Gibbs (full lines) and Boltzmann (dashed lines) Helmholtz free 
energies as a function of z (arbitrary units, N = 200). 

(a) Paramagnetic regime (T = 1.25 x T c ): z = 0 is the only stable equilibrium 
state, both for Boltzmann and for Gibbs. 

(b) Ferromagnetic regime (T = 0.75 x T c ): In addition to the two ferromagnetic 
states ±z we i SS , Gibbs predicts the permanence of a stable paramagnetic state at 
2 = 0 . 
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In conclusion, the difference between Boltzmann and Gibbs in a finite-size ferro- 
magnet, at zero external field, is that the former predicts an unstable paramagnetic 
state z = 0 in the ferromagnetic phase (T < T c ), while the latter predicts the persis¬ 
tence of a stable paramagnetic phase, cohexisting with the two ferromagnetic states 
±z we i ss - In other words, Boltzmann predicts what is called a ’’spontaneous symme¬ 
try breaking” at T = T c , in which an arbitrary small external field makes the system 
fall in one of the two ferromagnetic states ±z we i SS , at T < Tm- Gibbs, in contrast, 
predicts a lower limit B m of the external field, below which the paramagnetic state 
persists, even in the ferromagnetic phase. The calculation of B m = 2Jp m /po fol¬ 
lows from eq.n (20), by looking for the small extremants in the presence of a small 
field: 


d^cjz, p ) 
dz 


= 0 


2(z + p) = sign(z) ZG for \z\ < \ —jr: . 


For this equation to have two solutions, for both signs of z (which preserves the 
minimum in between), it is necessary that 

I . _ ZG _ T 

p Pm 2 T c V2Nn ' 

This shows that in a finite-size magnet, realized by elementary units of magnetic 
moment po, the symmetry breaking predicted by Gibbs is not ’’spontaneous”, but 
involves a lower limiting value of the external magnetic field: 


\B\ > B m = 


T 


PoV2Nt> 


{T < T m ) - 


(23) 


vanishing with the inverse square root of N. The decrease of B m (N) with increas¬ 
ing N corresponds to the ’’smoothing out” of the cuspid minimum. The 1/y/N- 
discrepancy between Boltzmann ( B m i n = 0) and Gibbs [B m in oc 1/y/N) is manifest 


from eq.n (23). 


The advantage of using a ferromagnetic material for testing the validity of Gibbs’ 
entropy is that the recent developments of micro and nano-physics make it possi¬ 
ble to prepare magnetic particles, containing controllable (and small) numbers of 
magnetic moments. Actually, size-effect measurements on such micro-magnets have 
become current in recent years |l5| . In the present case, one could measure B m 
with high sensistivity magnetometers, for different particles’ size. If SQUID magne¬ 
tometers were to be used, at the largest operative temperature of about 10 K, with 
Po ~ PBohr ) equation ( [23] ) yields B m i n « 10~ 6 G for a cube millimeter of magnetic 
sample, with a density N/V ~ 10 23 cm -3 typical of metals. 


5 Conclusions 

The contrast between Boltzmann’s and Gibbs’ approach to statistical thermody¬ 
namics is a long standing question, that comes to the light, time to time, since 
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about one century, in different contexts, but with the same underlying leitmotiv: 
which states are available to the system, and the way they must be counted. All 
this obviously reflects on the entropy. Dunkel and Hilbert, in ref. [T (see also 
ref. [132]), add a new element of discussion, by refreshing an almost forgotten defini¬ 
tion of entropy, due to Gibbs (eq.n ([3])), that forbids negative absolute temperature 
(NAT), in contrast to the current expression ([!]), usually attributed to Boltzmann. 
The arguments in ref. [I| have raised a discussion on the validity of eq.n ([3]), as 
an alternative to eq.n 0 [2109 m 0 EE [7] . The two expressions differ in the phase 
space regions uniformly occupied by the system: a surface of constant energy E, 
for Boltzmann; a volume containing all states with energy less than E, for Gibbs. 
While Boltzmann’s definition cannot apply but to an isolated, micro-canonical sys¬ 
tem, Gibbs entropy seems more appropriate for canonical systems, whose energy 
can fluctuate, due to the heat exchanges with a thermal bath. This, however, con¬ 
trasts with what is claimed in ref.s Hi nusnn]. Hence, in addition to the question 
about which of the two expressions is correct, one should also explore the possibil¬ 
ity that both expressions are correct, but refer, respectively, to micro-canonical and 
canonical constraints (as stressed in ref.s mil in ed , or vice versa. 

Since Gibbs and Boltzmann entropies yield the same results in the thermo¬ 
dynamical limit N —> oo, any possible measurable consequence of what precedes 
results in a size effect and is thereby far from easy to exploit. A superficial read¬ 
ing of eq.n (14) in ref. T] might lead to the conclusion that the smallness of the 
differences between eq.n 0 and (0 vanish as 1/N (the inverse of an extensive 
quantity like the heat capacity). If so, any experimental test, such as the ’’minimal 
quantum thermometer” suggested in ref. pQ, would become extremely difficult, if 
not impossible. 

The aim of the present work was to explore the possibility of concrete exper¬ 
imental procedures, deciding the winner (if any) of what we playfully called the 
Boltzmann vs Gibbs match, or claiming that the two opponents actually play on 
different playgrounds, if Boltzmann and Gibbs entropies would result to refer to 
different constraints. 

The study of two-level systems, both in the form of gases (Section 0, and of 
interacting Ising spins (Sections 0 and Appendix B), shows that the smallness fac¬ 
tor to be accounted for is not 1/iV, but 1 /VN, which makes any experimental 
procedure much more feasible. In particular, it has been seen that some cube cen¬ 
timeters of 3 He, with density about 10 23 cm -3 , at a temperature of tens Kelvins, 
under the action of a magnetic field of 10 - 'G, could be an appropriate system for 
an experimental test, possibly revealing the most striking effect of Gibbs entropy, 
i.e. a thermally activated overpopulation of the upper level. As a possible alterna¬ 
tive, it is found that the effect supporting Gibbs’ entropy in a Weiss ferromagnet 
(Ising model) should be a divergent heat capacity at a critical temperature Tm 
(lower than the standard Curie temperature), in contrast to the finite discontinu¬ 
ity predicted by Boltzmann’s picture at the Curie temperature T c (Figure 3). In 
a canonical process (Section 0, Gibbs’ entropy would be also responsible for an 
anomalous persistence of the paramagnetic regime, below the critical temperature 
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Tm, that should be destabilized by an external magnetic field proportional to l/VN 
(eq.n (231). This should be observed in canonical sub-millimetric metallic parti¬ 
cles, with controllable number of interacting magnetic moments, under the action 
of weak external magnetic fields (« 10 -6 G). 

The orders of magnitude involved in the measurements sketched above look 
fairly accessible to concrete experiments, which opens the possibility of an impor¬ 
tant advance in the thermodynamics of small systems, and in the genuine founda¬ 
tions of statistical thermodynamics. 


Post-amble 

As anticipated in the preamble, now I try to explain the reasons that lead one 
of the referees to reject the publication in PRE. 

The referee denies that equation ([2]), which is the cornerstone of all the question, 
is actually equivalent to the minimization of the Helmholtz free energy ’I' with 
respect to, say, the magnetization (or any other state parameter), at fixed V and 
T. The referee claims that, while this is true for Boltzmann’s entropy, it might be 
false for Gibbs’ entropy (giving no counter-example, however). In the lack of a first 
principle statistical derivation, says the referee, the minimization of d/ is an ad hoc 
procedure, and any result supporting it might be a pure coincidence. 

As shown by what precedes, my point is that equation © is manifestly equiv¬ 
alent to find the extremants of T, which is nothing but a mathematical outcome. 
On using just equation ([2]) (and nothing else), I show that, in the case of Weiss fer¬ 
romagnetism, some of those extremants behave as true thermal equilibrium values, 
while others behave like non equilibrium values (see Fig.2), yielding, for example, a 
negative heat capacity. Of course, the ‘good’ extremants and the ‘bad’ ones are the 
minima and the maxima of T, respectively (see Fig.s 5). As far as I can see, the 
only problem with this picture is the opinion , shared by many Gibbs’ supporters, 
that Gibbs’ entropy should apply to microcanonical systems only. If so, the notion 
of Helmholtz free energy is questionable, since it is well known that minimizing 
>F is a canonical equilibrium procedure. This is aknowledged as a possible pend¬ 
ing question to be explored, hopefully experimentally, which is right the aim of the 
present work: suggesting experimental tests to check the validity of two contrasting 
definitions of the entropy. 

Paradoxically, the efforts made in the (vain) attempt to remove the referee’s 
objections turned into some advantages. Thank to them, indeed, the work has 
been considerably improved. For example, the divergency of Cq at Tm (Fig. 3), 
that is an important physical aspect, had escaped to my attention in the first 
version. 
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A Appendix 

For TV 1, the contributes to the integral 


1 = 


J g-T[T>( z )+ a (®- 2 )+ fc (a:-2 : ) 2 + '"] c J a 


in eq.n © come from the minimum of (j)[x) (positive and symmetric) in the inte¬ 
gration interval [—1, z(E)\. So, cf>(x ) can be approximated by the first three terms 
in eq.n (5c), which yields: 


with 


rO 


/ e -f(ay+by 2 ) d y = J_ , 

/— (1+2:) -^AT 7- 


rO 


l n = 



(1 ~\~z)Ln 


N 


A " = a 'i 2b - 


-(Q 2 +A„Q) d g ; ( A . la ) 


(A.lb) 


If E < 0, one has z(E) < 0 in turn (eq.n ( [5b] )) and An < 0. If |Ajv| is large, Q 2 
can be neglected in the exponent of the integrand in eq.n (A.la), which yields: 


e f ~ _L f 
Ln J -o 


o — A/> 


°d Q = 


i(z)\N 


for |Ajv| > 1 . 


(A.2) 


On replacing expression (A.2) in eq.n (|6| ), one gets eq.n (7a), with the same validity 
condition (recall eq.ns (5c) and ( A.lb| 


If E > 0 (z, An > 0), the integral in eq.n (A.laI can be calculated by completing 
the square in the exponent: 


.£*(*)/ ; 


A%/4 /-Aiv/2 


e N 


e n 


/4 v^ 


Ln J ~ 


oo e ~ RdR= 2 Ln 


1 + erf 


An 

2 


(A.3) 


If, in particular, one takes 2 <C 1, equations (A.3), (5c) and (A.lb) lead to the 
expression in eq.n (7c). 

The integral 


J-T \4>{z)+a(x-z)+b{x-z ) 2 +---] da 


>\*\ 


(A.4) 


in eq.n © can be calculated with the same method, for |Ajy| 1 (eq.n ( |21a| ) ) and 
|Ajv| < 1 (eq.n (|21b|) ). 
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